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Abstract. We derive here a variant of the 2D Green-Naghdi equations that 
model the propagation of two-directional, nonlinear dispersive waves in shallow 
water. This new model has the same accuracy as the standard 2D Green- 
Naghdi equations. Its mathematical interest is that it allows a control of the 
rotational part of the (vertically averaged) horizontal velocity, which is not 
the case for the usual Green-Naghdi equations. Using this property, we show 
that the solution of these new equations can be constructed by a standard 
Picard iterative scheme so that there is no loss of regularity of the solution 
with respect to the initial condition. Finally, we prove that the new Green- 
Naghdi equations conserve the almost irrotationality of the vertically averaged 
horizontal component of the velocity. 



1. Introduction 

1.1. General setting. The water-waves problem, consists in studying the motion 
of the free surface and the evolution of the velocity field of a layer of fluid under the 
following assumptions: the fluid is ideal, incompressible, irrotationnal, and under 
the only influence of gravity. Many works have set a good theoretical background for 
this problem. Its local well-posedness has been discussed among others by Nalimov 
[T5] , Yosihara [T5] , Craig [5] , Wu jT5J US] and Lannes [5] ■ Since we are interested 
here in the asymptotic behavior of the solutions, it is convenient to work with a non- 
dimcnsionalizcd version of the equations. In this framework, (see for instance [2] 
for details), the free surface is parametrized by z = e£(t, X) (with X = (x, y) £ K 2 ) 
and the bottom by z = — 1 + @b(X). Here, e and f3 are dimcnsionlcss parameters 
defined as 



where a sur f is the typical amplitude of the waves and 6(, ott is the typical amplitude 
of the bottom deformations, while h is the depth. One can use the incompress- 
ibility and irrotationality conditions to write the non-dimensionalized water-waves 
equations under Bernoulli's formulation, in terms of a velocity potential ip associ- 
ated to the flow, where ip(t, .) is defined on fit = {(X, z),—l+/3b(x) < z < e((t, X)} 
(i.e. the velocity field is given by U = Vx,z<p) ■ 
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(1) 



H%tp + fxd 2 ip + d 2 z <p = 0, 
d z tp - ■ V(p = 

9tC-i(M£VC-V^ + ^) = 0, 
dw+ l -{e\V v \ 2 + ^{d z v?) + t = Q 



at -1 + /?&< 2 < e£, 

at z = -l+fib, 

at z = , 

at z = e(, 



where V = (d x , d y 



(di,d2) T - The shallowness parameter jj, appearing in this 



set of equations is defined as 



A 2 ' 



where, A is the typical wave-length of the waves. Making assumptions on the size 
of e, /3, and \i one is led to derive (simpler) asymptotic models from (JTJ) . In the 
shallow- water scaling (// <C 1), one can derive (when no smallness assumption is 
made on e and f3 ) the standard Green-Naghdi equations (see £12.11 below for a 
derivation and [2] for a rigorous justification). For two-dimensional surfaces and 
over uneven bottoms these equations couple the free surface elevation ( to the 
vertically averaged horizontal component of the velocity, 



(2) 



v(t,X) = 



1 



1 + e( - f3b 
and can be written as: 

d t ( + V • (hv) = 0, 
(h + fj,T[h, /36]) d t v + hV£ + 



l+/3b 



Vip(t,X,z)dz, 



(3) 



< 



+ fie{-\7[(h 3 (d 1 v-d 2 v J 



where h = 1 + e£ - 
(4) T[h,sb]W 



fib, v = (V u V 2 ) T 
1 



■e(h + nT[h,l3b])(v V)« 
+ (V-«) 2 )]+5R[/i,/36](«)} 

- = (-^ 2 ,Vi)- L ,and 

/3 r 



V(/i 3 V • VF) + ^[V(A a V6 ■ WO - /i 2 V6V ■ W] 
+f3 2 hVbVb ■ W, 
while the purely topographic term j3b](v) is defined as: 

Pi 



(5) 



R[h,)86](t;) = ^V(h 2 (V 1 2 d 2 1 b + 2V 1 V 2 d 1 d 2 b + V 2 2 dlb)) 

+f3h 2 (d 1 v ■ d^ 1 - + (V • v) 2 )Vb 
+(3 2 h(V 2 d' 2 b + 2V 1 V 2 d 1 d 2 b + V 2 d 2 2 b)Vb. 



A rigorous justification of the standard GN model was given by Li QT] in ID 
and for flat bottoms, and by B. Alvarez-Samaniego and D. Lannes [2] in the general 
case. This latter reference relies on well-poscdness results for these equations given 
in [3] and based on general well-posedness results for evolution equations using a 
Nash-Moser scheme. The result of [3] covers both the case of ID and 2D surfaces, 
and allows for non flat bottoms. The reason why a Nash-Moser scheme is used there 
is because the estimates on the linearized equations exhibit losses of derivatives. 
However, in the ID case, such losses do not occur and it is possible to construct a 
solution with a standard Picard iterative scheme as in PH E] with flat and non flat 
bottoms respectively. But, this is not the case in 2D, since for instance, the term 
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d\v ■ 82V 1 - is not controled by the energy norm | • \y naturally associated to ([3]) 
(see [3]), 

|(C,«)|^ = |CI^ + I«I^+m|V-«| 2 h ,. 

This is the motivation for the present derivation of a new variant of the 2D 
Green-Naghdi equations (|3|). This variant has the same accuracy as the standard 
2D Green-Naghdi equations (J3j> (sec i 32.2l bclow for a derivation) and can be written 
under the form: 

- d t c + v • (hv) = o, 

[h + fi(T[h, Pb] - V^curl )) d t v + hV( 

+ e(h + n(T[h, /3b] - V^curl )) (v ■ V)v 

+ ^e^V[h 3 (d lV ■ 82V 1 - + (V • v) 2 )} + K[h,pb](v)} = 0, 

where V -1 = (— d y ,d x ) T , curlw = d\V2 — 82V!, while the linear operators T[h,eb] 
and 3?[/i, (3b] are defined in (jU) and (0. 

The reason why the new terms (involving V^curl ) do not affect the precision of 
the model is because the solutions to © are nearly irrotational (in the sense that 
curlw is small). This property is of course satisfied also by our new model ([6]). 
The presence of these new terms allows the definition of a new energy norm that 
controls also the rotational part of v. Consequently, we show that it is possible to 
use a standard Picard iterative scheme to prove the well-poscdncss of (j6|), so that 
there is no loss of regularity of the solution with respect to the initial condition. 

1.2. Organization of the paper. We first recall the derivation of the standard 
2D Green-Naghdi equations in Section |2~T1 while in Section l2~2l we derive the new 
model ([6]) . We give some preliminary results in Section I3.1( the main theorem 
which proves the well-posedness of this new Green-Naghdi system is then stated 
in Section [3721 and proved in Section [3731 Finally, in Section [3741 we prove that J6|) 
conserves the almost irrotationality of v. 

1.3. Notation. We denote by C(Ai, A2, ...) a constant depending on the parameters 
Ai, A2, ... and whose dependence on the Xj is always assumed to be nondecreasing. 
The notation a < b means that a < Cb, for some nonegative constant C whose 
exact expression is of no importance (in particular, it is independent of the small 
parameters involved). 

Let p be any constant with 1 < p < 00 and denote L p = L P (BL 2 ) the space of all 
Lebesguc-measurablc functions / with the standard norm 

|/|l* = ( f \f(X)\ p dX) 1/p <oo. 

JR 2 

When p = 2, we denote the norm \ ■ \l 2 simply by | ■ I2. The inner product of any 
functions fi and f2 in the Hilbert space L 2 (M. 2 ) is denoted by 

(h,h)= j Ji{X)f 2 {X)dX. 

The space L°° = L°°(1R 2 ) consists of all essentially bounded, Lebesgue- measurable 
functions / with the norm 

|/|l°o = esssup|/(X)| < 00. 
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We denote by W^°° = W 1 ' 00 ^ 2 ) = {/ e L°°,Vf G (L°°) 2 } endowed with its 
canonical norm. 

For any real constant s, H a = H S (M. 2 ) denotes the Sobolev space of all tempered dis- 
tributions / with the norm = |A S /|2 < oo, where A is the pseudo-differential 
operator A = (f - A) 1/2 . 

For any functions u = u(x,t) and v(x,t) defined on R 2 x [0,T) with T > 0, we 
denote the inner product, the L p -norm and especially the L 2 -norm, as well as the 
Sobolev norm, with respect to the spatial variable X, by (u,v) = (u(-,t),v(-,t)), 
\u\lp = \u(-,t)\ L p, \u\ L 2 = \u(-,t)\ L 2 , and \u\ H * = \u(-,t)\ H *, respectively. 
Let C fc (R 2 ) denote the space of fc-times continuously differentiable functions and 
C^°(R 2 ) denote the space of infinitely differentiable functions, with compact sup- 
port in R 2 ; we also denote by C£°(R 2 ) the space of infinitely differentiable functions 
that are bounded together with all their derivatives. 

For any closed operator T defined on a Banach space X of functions, the commu- 
tator [T, f] is defined by [T, f]g = T(fg) - fT(g) with /, g and fg belonging to the 
domain of T. 

We denote v 1 - = (— V 2 , V\) T and curlw = d\Vi — d^Vx where v = (Vi,V2) T . 



2. Derivation of the new Green-Naghdi model 

This section is devoted to the derivation of a new Green-Naghdi asymptotic 
model for the water-waves equations in the shallow water (p -C 1) of the same 
accuracy as the standard 2D Green-Naghdi equations (0). 

2.1. Derivation of the standard Green-Naghdi equations ([3]). We recall 
here the main steps of jTD] for the derivation of the standard 2D Green-Naghdi 
equations ([3]). In order to reduce the model (|T|) into a model of two equations we 
introduce the trace of the velocity potential at the free surface, defined as 

and the Dirichlet-Neumann operator C/ M [eC,/36]- as 

with if solving the boundary value problem 

{ IJ^lf + IJ-dyf + d 2 z f = 0, 

(7) \ d n ip u _ 1+[lb =0, 

As remarked in [19l [6] , the equations ([l) are equivalent to a set of two equations 
on the free surface parametrization £ and the trace of the velocity potential at the 
surface ip = y|. = „ c involving the Dirichlet-Neumann operator. Namely 

(8) 3 , ! ^,v7/i2 (^mK,W + v( £ C)-v^) 2 

M + C + 2^ - ^ 2(1 + £ W) = °' 

It is a straightforward consequence of Green's identity that 

-s M K/3*# = -v-(M, 
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with h = 1 + e£ — (3b and v = 1+ J-_p b J^_ 1+ g b V(p(t,X,z)dz. Therefore, the first 
equation of © exactly coincides with the first equation of ([3]). In order to derive 
the evolution equation on v, the key point is to obtain an asymptotic expansion of 
Vip with respect to \i and in terms pf v and £. As in [10] . we look for an asymptotic 
expansion of ip under the form 

N 

(9) (p app ^^2fi J (pj. 

3=0 

Plugging this expression into the boundary value problem ([7]) one can cancel the 
residual up to the order 0(n N+1 ) provided that 

(10) Vj = 0, N, 8 2 zVj = - tftpj-x 

(with the convention that (p—% = 0), together with the boundary conditions 



(11) Vj = 0,..,jV, 



d z ifj = (3Vb • V^_i Ue= _ i+/js 



(where <5oj = 1 if J = and otherwise). 

By solving the ODE (fT0|) with the boundary conditions (fTTj) . one finds (see [10] ) 

(12) <p = V, 

(13) tpi = {z-e()(-^(z + e()-l + pb)AiP + f3{z-e()Vb.Vip. 



According to formulae (fT2j) . (|13[) . the horizontal component of the velocity in the 
fluid domain is given by 

V(z) = VMz) + V<pi{z) + 0(m 2 )- 

The averaged velocity is thus given by 

1 



h 



or equivalcntly 



(V^o(z) + V^iOb)) dz + 0(n 2 ), 
-1+Bb 



VV>-/^T[M&]VV> + 0(/i 2 ), 



and thus 

(14) Vip = v + ^T[h,f3b}Wt/j + 0(n 2 ), 

where T[h, /3b] is as in @. As in [TU], taking the gradient of the second equation 
of ([HI); replacing Vtp by its expression (jT4")) and —Gfi,[£(, (3b]tp by —V ■ (hv) in the 
resulting equation, gives the standard Green-Naghdi equations (after dropping the 
0(n 2 ) terms), 

f d t C + V • (hv) = 0, 

(h + fiT[h, (3b})d t v + hVQ + e(h + nT[h, (3b}) (v • V)« 

+ /ie{^V[/i 3 (9n; • d 2 v^ + (V • v) 2 )] + 3t[h, /3b](v)} = 0, 

where h = 1 + e( — (3b, v = (Vi, V2) T , v 1 - = (— V2, Vi) T , and the linear operators 
T[h,eb] and ^R[h,(3b] being defined in (gD and ©. 
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2.2. Derivation of the new Green-Naghdi system (|6|). It is quite common in 
the literature to work with variants of the Green-Naghdi equations ((6]) that differ 
only up to terms of order 0(/j 2 ) in order to improve the dispersive properties of the 
model (sec for instance [Mill]) or to change its mathematical properties [12] . Our 
approach here is in the same spirit since our goal is to derive a new model with 
better energy estimates. 

In order to obtain the new Green-Naghdi system ([6]), let us remark that, from the 
expression of v one has 

VV> = v + ^T\h,/3b}V^ + 0( f i 2 ); 
and since v = VV> + 0(/i), one gets the following Lemma: 

Lemma 1. Let v be the vertically averaged horizontal component of the velocity 
given above, one obtains 

curl(9 t v = O(p), 
(15) curl(VV)u = O(p). 

Remark 1. For the sake of simplicity, we denote by O(n) any family of functions 
(/m)o</i<i such that {j^f^) 0<tl<1 remains bounded in L°°([0, ~], H r (R 2 )), for some 
r large enough. 

Proof. By applying the operator (curl<9t)(-) to the identity v = S7ip + 0(/j,) one gets 

curlew = 0(/i). 

In order to prove the second identity of (|15|) . replace v = \7tp + 0([a) in (v ■ V)i> 
and apply the operator curl (•) to deduce 

curl (v ■ V)v = 0{n). 

□ 

Using Lemma [J the quantities /iV^curlc^w and fjiV curle(u • V)i> are of size 
0(/x 2 ), which is the precision of the GN equations We can thus include these 
new terms in the second equation of ([3]) to get 

' dtC + V • (hv) = 0, 

(h + n(T[h, (3b] - V^curl ))d t v + hV( 

+ e(h + n(T[h, f3b] - V^curl )) (v • V)w 

+ f Me{^V[h s (d 1 v ■ 32V 1 - + (V • v) 2 )] + St[h,pb](v)} = 0(/i 2 ), 

(the linear operators T[h,eb] and $t[h,j3b] being defined in ^ and (O). 

Remark 2. We added the quantity /xV^curldiV = 0(fi 2 ) to the standard GN 
equations to obtain an energy norm \ ■ \x a - 

l(C,«)&. = Kl 2 ff * + 14. +m|v • V\* H . + »\cuv\v\ 2 Ha . 

The last term is absent from the energy \ ■ |y« associated to the standard GN equa- 
tions (0) (see We will see in the next section that it allows a control of the 
term d\v ■ 82V 1 - . 
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Remark 3. The bilinear operator $t[h,/3b](v) only involves second order deriva- 
tives of v while third order derivatives of v have been factorized by h + /i(fT[h,eb} — 
V^cwrZ ) . The fact that the operator $t[h, f3b](v) does not involve third order deriva- 
tives is of great interest for the well-posedness of the new Green-Naghdi model. 

3. Mathematical analysis of the new Green-Naghdi model 

3.1. Preliminary results. For the sake of simplicity, we take here and throughout 
the rest of this paper (/? = e) and we write 

% = h + n{T[h,Eb]-V ± cm\). 

We always assume that the nonzero depth condition 

(16) 3 h mm > 0, inf h > h rnin , h = 1 + s(( — b) 

YGR 2 

is valid initially, which is a necessary condition for the new GN type system ([6]) to 
be physically valid. We shall demonstrate that the operator T plays an important 
role in the energy estimate and the local well-posedness of the GN type system ([6]) . 
Therefore, we give here some of its properties. 

The following lemma gives an important invertibility result on T and some prop- 
erties of the inverse operator 

Lemma 2. Let b e C b °°(R 2 ), t > 1 and C € H to+1 {R 2 ) be such that fH) is 
satisfied. Then, the operator I has a bounded inverse on (L 2 (IR 2 )) 2 , and 

(i) For all < s < to + 1, one has 

|T~7l^ + v^iv ■ srVk- + v^lcurlT- 1 /!^ < C(-^—, \h - l\ Ht0+1 )\f\ Hs , 
and 

^- l Vg\ Hs + v^lS-'V^lffa < C{-^—, \h - l\ Ht0+ ,)\g\ H s. 

(ii) Ifs>t + l and(e H S (R 2 ) then 

II £ _1 \\{H°) 2 ->(H°) 2 +VM II 'J _1 V \\(H>) 2 ^(H>) 2 +V / M II T _1 V \\(H') 2 ^{H<>)2< c s, 

and 

II v ' ^ 1 ll(H s ) 2 ^(ff s ) 2 +VM II curlT -1 j|(^) 2 ^(H«) 2 < c s , 
where c s is a constant depending on ■— , \h — 1\h*> and independent of (\i,e) 

e (o,i) 2 . 

Remark 4. Here and throughout the rest of this paper, and for the sake of sim- 
plicity, we do not try to give some optimal regularity assumption on the bottom 
parametrization b. This could easily be done, but is of no interest for our present 
purpose. Consequently, we omit to write the dependence on b of the different quan- 
tities that appear in the proof. 

Proof. It can be remarked that the operator T is L 2 self-adjoint; since, one has 
(h + n(T[h, eb] - V^curl )v,v\ = {hv, v) 

+^(/i(-^=V ■ v - e^-Vb ■ vj, -^=V • v - £-^V6 • vj + ^-{hVb ■ v,Vb-v) 
+/x(curlii, curlw), 
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and using the fact that infR2 h > h m i n , one deduces that 

(1v,v) > E{eb,v), 

with 

E(eb,v) := h min \v\ 2 2 + fih min \^=V ■ v - e^-Vb ■ v\% + M£ ^ mm [V6 ■ v\ 2 + n\cwlv\l 

proceeding exactly as in the proof of the Lemma 1 of [7] it follows that T has an 
inverse bounded on (L 2 (R 2 )) 2 . 

For the rest of the proof, One can proceeding as in the proof of Lemma 4.7 of [3J, 
to get the result. □ 

3.2. Linear analysis of ([6]). In order to rewrite the new GN type equations © 
in a condensed form, let us write U = (C,v T ) T , v = (Vi, V 2 ) T and 

Q = h 3 (-d 1 vd 2 v ± -(V-v) 2 ). 

We decompose now the 0(e/x) term of the second equation of ([6|) under the form 

-^7Q + R[h,£b]{v) = Rx[U]v + r 2 {U), 

with for all / = (F 1 ,F 2 ) T 

Ri[U]f = -^Q[U]f+^(h 2 (V 1 F 1 d 2 b + 2V 1 F 2 d 1 d 2 b + V 2 F 2 d 2 b)) 

-£h 2 (-d lV ■ d 2 f^ - (V • w)(V • f))Vb; 

r 2 (U) = e 2 h(V?d 2 b + 2V 1 V 2 d 1 d 2 b + v£d$b)Vb, 

where 

Q[U]f = h 3 (-d lV ■ 8 2 f^ - (V • «)(V • /)), 

(in particular, Q = Q[U]v). The new Green-Naghdi equations (|6|) can thus be 
written after applying T _1 to both sides of the second equation in (jSJ) as 

(17) d t U + A[U]U + B(U) = 0, 

with U = (C, V u V 2 ) T , v = (Vi, V 2 ) T and where 



(18) A[U] 
and 

(19) 



ev ■ V hV- 



This subsection is devoted to the proof of energy estimates for the following initial 
value problem around some reference state U_ = (£, Vi, V 2 ) T : 

(20) [ d t U + A[U]U + B(U) = 0- 

We define first the X s spaces, which are the energy spaces for this problem. 
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Definition 1. For all s > and T > 0, we denote by X s the vector space H S (M. 2 ) x 
(H S+1 (R 2 )) 2 endowed with the norm 

VU=((,v)e X s , \U\\ S := IClfr- + M(ff-) a + m|V ■ v\ 2 H . + v\cwlv\ 2 HS , 

while Xj, stands for C([0, -j]; X s ) endowed with its canonical norm. 

We define the matrix S as 

( 21 > s=(ll). 

with h = l + e((-b) and X = h + fi(T\h, eh] - V J -(VA )). A natural energy for the 
IVP ((201) is given by 

(22) E s {Uf = (A S U,SA S U). 

The link between E S (U) and the A s -norm is investigated in the following Lemma. 

Lemma 3. Let b E C fc °°(R 2 ), s > and ( e W 1 '°°(R 2 ). Under the condition Ho)) . 
E S (U) is uniformly equivalent to the | • \x°-norm with respect to (//, e) € (0, l) 2 : 

E s (U)<C(\h\ wl ,oo)\U\ x *, 

and 

\U\ X s<C(-^-)E%U). 

nrain 

Proof. Notice first that 

E S (U) 2 = |A s Cll + (A s v,TA s w), 

one gets the first estimate using the explicit expression of °\£, integration by parts 
and Cauchy-Schwarz inequality. 

The other inequality can be proved by using that mi xe ^2 h > h m in > and pro- 
ceeding as in the proof of Lemma [2j □ 

We prove now the energy estimates in the following proposition. It is worth 
insisting on the fact that these estimates are uniform with respect to e, [i € (0, 1); 
since the control of the s + 1 order derivatives by the X s -norm disappears as /i — > 0, 
the uniformity with respect to /i requires particular care (see the control of B^q in 
the proof for instance), but is very important for the application since fi <C 1. 

Proposition 1. Let b e C£°(R 2 ), t > \, s > t + 1. Let also U = (C, Vl, V2_) t 
£ Xj, be such that dtU_ £ X^T and satisfying the condition S16\) on [0, Then 
for all Uq G X s there exists a unique solution U = {(,Vi,V 2 ) T € Xf, to and 
for all0<t< 

E s (U(t)) <e EXTt E s (U )+e f e 6 ^^^ C{E S {U_){t'))dt' , 

Jo 

for some X T = A T (sup < t < T/e E s (U(t)), sup < t < T/e \d t h{t)\ L ^) . 

Proof. Existence and uniqueness of a solution to the IVP (|20)) is achieved by using 
classical methods as in appendix A of [7] for the standard ID Green-Naghdi equa- 
tions and we thus focus our attention on the proof of the energy estimate. For any 
A 6 R, we compute 

e eXt d t {e- £Xt E a {U) 2 ) = -e\E s (U) 2 +d t {E s (U) 2 ). 
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Since 

E S {U) 2 = (A S [/,SA S £7), 
and U = (C, Vi, V 2 ) T , v = (V U V 2 ) T , we have 

(23) d t (E s (Uf) = 2(A s (,A s ( t ) + 2(A s v,SA s v t ) + (A s v,[d t ,%A s v). 
One gets using the equations (|20j) and integrating by parts, 

^e eM d t {e- eXt E s {U) 2 ) = -y£ s (f/) 2 - (SA[U]A S U, A S U) 

(24) -([A S ,A[U]]U,SA S U) - (A S B(U), SA S U) + i(A s v, [<9 t , T]A s i;). 

We now turn to bound from above the different components of the r.h.s of (|2~H) . 
• Estimate of (SA[C/]A S C7, A S U). Remarking that 

SA ^ = ( ~hV %{ev ■ V)+ enRi [U] ) ' 

we get 

{SA[U]A S U,A S U) = (ev-VA s CA s () + (hV -A s v,A s () 

+(WA s (,A s t.) + ((%(ev • V) +£^i?i[C/])A s v,A s w) 
=: Ai + A 2 + A 3 + A A . 

We now focus on the control of (A,-)i<j<4. 

— Control of A\. Integrating by parts, one obtains 

Ax = (ev ■ VA S C, A S C) = ~(eV • «A S C, A S C), 
and one can conclude by Cauchy-Schwarz inequality that 

|Ai| <eC(\V -v\ L ~)E s (U) 2 . 

— Control of A 2 + A3. First remark that 

\A 2 + A 3 \ = \(Vh-A s v,A%)\ < \Vh\ {L ~ r E s {U) 2 - 

we get, 

\A 2 +A 3 \ <eC(\Vh\ {L ^E s (U) 2 . 

— Control of A4. One computes, 

Ai = e{%{(v-^)A s v),A s v)+{e^R 1 [U]A s v 1 A s v) 
=: A 41 +A 42 . 

Note first that 

A u = e(h(v-V)A s v,A s v) + ^(h 3 V-(v-V)A s v),A s V-v) 

O 

2 2 

-^-y(h 2 Vb-{v- V)A s v, A S V ■ v) - ^(h 2 VbV ■ (v ■ V)A s v, A s v) 
+e 3 fj,(hS7bS7b- (v-V)A s v,A s v) + e/i(curl (v ■ V)A s v), A s curlw); 
remark also that 

(curl (v ■ V)A s w), A s curlu) = ~ ((curl A s v, d^curl A s v) + (V A A s v, d 2 V 2 y A A s v) 

+(curlA s w,9n) • VA s y 2 ) - (curl A s v, d 2 v • VA s lq), 
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and that, for all F and G smooth enough, one has 

((G • V)v, F) = -(«, F V • G) - (v, (G ■ V)F). 
By using successively the above identities, the following relation (f2"5j) 

(25) iVFrll + |VF 2 |2 = |V • F\\ + \cuv\F\ 2 , 

integration by parts and the Cauchy-Schwarz inequality, one obtains directly: 
\A ll \<eC(\^ (W i,~ ) *,\C\ w i.~)E'(U) 2 . 

For A42, remark that 

|^ 42 | = |(^i[C/]A s «,A s «)| 

- \en{Q[U]A 8 v, A S V • v) - tie 2 {^{-dyv ■ d 2 A s v ± - (V • u)(V • A s v))Vb, A s i 

2 

-^{tk {]LiA s Vxdlb + 2V 1 A s V 2 did 2 b + V 2 A s V 2 dib), A S V • v) , 

where 

QU]f = ^(-div ■ d 2 f x - (v • «)(v • /)). 

we deduce that 

\A i2 \ < eC(\v\ {w i,^,\C\w^)E s (U) 2 . 
The estimates proved in A41 and A 42 show that 

\M <eC{\v\ (w i,~y,\Q w i,~)E a {U) 2 . 

• Estimate of ( [A s , A[U}] U, SA S U). Remark first that 

([A s , A[U]}U,SA S U) = ([A s , e «] • VC, A S C) + ([A s ,h]V- v,A s C) 

+ ([A S ,T _1 h}V(,%_A s v) + {[A s ,e(v ■ V)]v, TA s w) + e^([A s ,T^ 1 i?i[C/]]w,TA s i;) 

=: B 1 +B 2 +B 3 + B i + B B . 

- Control oiB l +B 2 = ([A s ,ev] ■ VC,A S C) + ([A s ,ftjV • v, A s (). 

Since s > to + 1, we can use the following commutator estimate (|26p (see e.g [9]) 

(26) |[A',JlG| a <|VF| H .-i|G|B.-i. 
to get 

\Bt+B 2 \ < eC{E s (U))E s {U) 2 . 

— Control of B4 = ([A s ,e(v ■ V)]u, TA s w). By using the explicit expression of T we 
get 

B A = ([A s ,ev ■ V]v,hA s v) + ^(V • [A s ,e(i;- V)]v,h 3 A S V ■ «) 

-y ([A s ,e(z; • V)]v,/i 2 V6A S V • «) + y([A s , e(v ■ V)]v, V(h 2 Vb-A s v)) 

+£ 2 ^([A S , e(v • V)]u, h V6V6 • A s v) + //(curl [A s , e(u • V)]u, A s curlu) 
:= B41 + B42 + B43 + -S44 + -B45 + -046. 
One obtains as for the control of B\ and i?2 above that 

\B 4j \<eC(E s (lLj)E s (U) 2 , j £{1,3,4,5}. 



The control of -B42 and B^& is more delicate because of the dependence on /i (recall 
that the energy E S {U) controls s + 1 derivatives of v, but with a small coefficient 
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^/JL in front of the derivatives of order 0{yfp)). For B 46 , we thus proceed as follows: 
we first write 

B 46 = ^(curl [A s ,e(v ■ V)]u, A s curlw) 

= m(9i[A s ,£_^9i]1/ 2 , A s curlv) + M (ai[A s ,e_^a 2 ]y 2 ,A s curl«) 
-fi(d 2 [A s ,eVidi}Vi, A s curlt;) - n{d 2 [A s , e\^d 2 ]Vi, A s curlv) 

:= £>461 + B 46 2 + -B463 + #464- 

Remarking that for all j S {1, 2} we have 

d j [A s ,f]g=[A s ,d j f]g+[A s ,f]d j g, 

and 

[A s ,fd j ]g = [A s ,f]d j g, 
we can rewrite B i6 i under the form 

#461 = K&IAVeH&lV&.A'curltO 

= ^([A s ,ed 1 V L ]d 1 V 2 ,A s cui\v) + /i([A s , eVi\dfV 2 , A'curlu). 

It is then easy to use the commutator estimate (|26[) in order to obtain 

|£? 4 <y| < sC(E s (U))E s (U) 2 , j € {1,2,3,4}. 

Similarly, B 42 is controled by eC(E s (V_))E s (U) 2 . The estimates proved in 
(-B4i)je{i,2,3,4,5,6} show that 

I-B4I < eC(E s (U))E s (U) 2 . 

- Control of B 3 = ([A S ,T _1 h]V(,%A s v). Remark first that 

£[A S ,S _1 ]WC - SA 5 ,^ 1 h]V( - [A s ,h]V(; 
moreover, since [A S ,T _1 ] = — %~ [A S ,T]£~ , one gets 

S[A S ,S _1 = -[A S ,2]S _1 WC + [A s , h]VC, 
and one can check by using the explicit expression of % that 

o 



^-V[A s ,/i 2 V6] • %- 1 hV(+ ^-[A s ,h 2 Vb]V • (T _1 ftVC) 



Integrating by parts yields 
S 3 = OAM^/ijVCAV) 



3 

+^ ( [A S , h 2 V6] • r'wc, A*V • t') + y 
-e 2 fi([A s , W6V& T ]2 _1 WC, A s «) + ([A s ,/i]VC,A s «) 



([A s ,Zi]T~ 1 WC, A s v) - |({[A s ,/i 3 ]V ■ (I _1 WC)}, A S V ■ «) 
^([A s ,/i 2 V6] -I-^VC.A'V -v) + y([A s ,/i 2 Vfo]V • (1 _1 W0,A S «) 
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One deduces directly from Lemma [21 the commutator estimate ([26]) . and Cauchy- 
Schwarz inequality that 

|S 3 | < C(-^—,\h-l\ H s) {(\Vh\ H .-i +e 2 v\h 2 VbVb T \ H s +^\hVb\ H -)\hV<;\ H . 

+(^|v/iV-i + ^|a 2 v6|h.)|wcu.-i + \vh\ H .- 1 \v<;\ H .-i)}\v\x: 

Finally, we deduce 

\B 3 \<eC(E s (U))E s (U) 2 . 

- Control of B 5 = s^([A s R^U^v^A^) . Let us first write 

%[^ s ^ Ri[US\v = -[A s ,Z\%- 1 R 1 [U]v + [A a ,Ri\U\]v 

so, that 



3 

_^ V {[A s ,^ 2 V6]-T- lD 
-e 2 /j,[A s ,hVbVb T ]'Z- 1 R 1 \U}v + [A s ,Ri[U]]v 



%[A S R\[U}\v = -[A s ,h]T 1 R 1 \U]v + ^{[A s ,h 3 ]V ■ 

2 V{[A s ,/i 2 V6] ■T- 1 -Ri[C/]u} + y[A s ,/i 2 V&]V • (^Ril^v) 



To control the term ( [A s , R\ [U]\ v, A s v) we use the explicit expression of Ri[U]: 
-etfi-diV ■ 8 2 f^ - (V • v)(V ■ f))Vb, 



Ri[U]f = -r^QMf + l^{h\V 1 F l dib + 2V 1 F 2 d 1 d 2 b + V2F 2 d I 2 b)) 



where, 

Q[U]f = h a (-a 1 vd a f ± - (V • w)(V • /)). 

As for the control of the term (V{[A S , h 3 ]\7 ■ Ri[U]v)} , A s v) we use the explicit 
expression of R\ [U_] , the relation (|25[) , the commutator estimate (|26|) and Lemma 
[2 Indeed, 

(V{[A s ,/i 3 ]V ■ (S-^itt/ju)}, A\-) = -^([A s ,/i 3 ]V ■ (T-^QIC/]^, A S V • w) 
| ([A s , Zi 3 ]V • T~ 1 V(fe a Q5yia?& + 2tfy 2 0ift6 + V^Vrfib)), A S V • «) 
- e([A s , h 3 ]V ■ ^ x (/i 2 (-5iu • d 2 A s v ± - (V ■ u)(V ■ w))V&), A S V ■ v). 
and thus, after remarking that 

\V > (g- 1 VQ\U]v)\ H .-i < I^VQItfH*., 
we can proceed as for the control of B$ to get 

|S 5 | <eC(E s (U))E s (U) 2 . 



Estimate of (A S B(U), SA S U). Note first that 

B(U) = 



e\7b ■ v 

-l 



e^~r 2 {U) 
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so that 

(A S B(U),SA S U) = {A^eVb-v),A s () + (A s ('Z- 1 r 2 {U)),'ZA s v) 

= (A s (eV6 • v), A S Q - £/i([A s ,T]T"V 2 (^), A s v) 
+e/i(AV 2 ([/),A s t>). 
Using again here the explicit expressions of %, r 2 (U_) and Lemma [H we get 

(A S B(U),SA S U) < eC(E s (U))E s (U). 
• Estimate of (A s v, [d t ,%]A s v). We have that 

(A s v,[d t ,%]A s v) = (A s v,d t hA s v) + ^(A S V •v,d t h 3 A a V • v) 

-^-(A s v,d t h 2 \7bA s \7 ■ v) - y (A S V • v,d t h 2 \7b ■ A s v) 
+e 2 n{A s v, d t hVbVb T A s v). 

Controlling these terms by sC(E s (U), \dth\L<=c )E S (U) 2 follows directly from a Cauchy- 
Schwarz inequality and an integration by parts. 

Gathering the informations provided by the above estimates and using the fact 
that the embedding H S (R 2 ) C W 1,00 (M. 2 ) is continuous, we get 

e eXt d t (e- eXt E s (Uf) < e(C(E s (U), \d t h\ L oo) - \)E S (U) 2 + eC(E s (U))E s (U). 

Taking A = At large enough (how large depending on sup tg j T] C(E s (U_), \dth\L°°) 
to have the first term of the right hand side negative for alH € [0, ~], one deduces 

V< e [0, -1, e eM d f {e- £Xt E s (U) 2 ) < eC(E s (U))E s (U). 

e 

Integrating this differential inequality yields therefore 

Vt e [0, -], E S {U) < e eXTt E s (U ) + e [ C{E S (U)(t'))dt' , 

£ Jo 

which is the desired result. □ 

3.3. Main result. In this subsection we prove the main result of this paper, which 
shows the well-posedncss of the new Green-Naghdi equations ([6]) for times of order 

Theorem 1. Let b g C & °°(R 2 ), i > I, s > i + 1. Let also U a = (Co,Vo) T e x * 
be such that ilo}) is satisfied. Then there exists a maximal T max > 0, uniformly 
bounded from below with respect to s,fj, £ (0,1), such that the new Green-Naghdi 
equations admit a unique solution U = (£,v T ) T G A" T with the initial 
condition {Coi v o) T an< ^ preserving the nonvanishing depth condition U6\) for any 
t 6 [0, )■ Ln particular if T max < oo one has 

|tf(t,.)|*.— >oo as t ^^-, 

£ 

or 

inf h(t,-) = inf 1 + e(((t,-) - &(•)) — >0 as t — > Tl0f±, 
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Remark 5. For 2D surface waves, non flat bottoms, B. A. Samaniego and D. 
Lannes [3] proved a well-posedness result for the standard 2D Green-Naghdi equa- 
tions using a Nash-Moser scheme. Our result only uses a standard Picard iterative 
and there is therefore no loss of regularity of the solution of the new 2D Green- 
Naghdi equations with respect to the initial condition. 

Remark 6. No smallness assumption on e nor /i is required in the theorem. The 
fact that T max is uniformly bounded from below with respect to these parameters 
allows us to say that if some smallness assumption is made on e, then the existence 
time becomes larger, namely of order 0(l/e). 

Proof. Using the energy estimate of Proposition [I] one proves the result following 
the same lines as in the proof of Theorem 1 in [7] , which is itself an adaptation of 
the standard proof of well-posedness of hyperbolic systems (e.g [HUT])- '-' 

3.4. Conservation of the almost irrotationality of v. To obtain the new 2D 
Green-Naghdi model ((6|) we used the fact that curlw = O(fi), we prove in the 
following Theorem that the new model ([6]) conserves of course this property. 

Theorem 2. Let t > 1, s > t + I, U = ((o,v?) T € X s+A with |curlu | ff = < 
LlC{\Uq\x-)- Then, the solution U = ((,v T ) T € X^+ 4 ax of the new Green-Naghdi 
equations {6]] with the initial condition (Co,Vq) t satisfies 

VO < T < T max , |curlv| L3O([0 i:] iff3) < fiC(T, \U \ X *+*)- 

Proof. Applying the operator curl (•) to the second equation of the model (|6|) after 

multiplying it by — yields 
h 

dtw + eV • (vw) — sfifi + g, 

where 

f l = -curl ( r T[h, eb] - -J-V^curl )v ■ Vw 
y h h ' 

-curl I j^V[/i 3 (div • + (V • v) 2 )} + R[h,e6](«)| 

and 

g = — itcurl (— T[h, eb] — — V^curl )d t v. 
h h 

From the identity curl(4W) = — £t^V" L (C — b) ■ W + ^curl W, we deduce that g 
can be put under the form 

9 = £Vfi + MCui'l(T V± dt w )- 
h 

with 

f2 = -^V ± ((-b)-T[h,sb]d t v. 
We have thus shown that w solves 

(j - jucurl (-V 1 )^Jd t w + eV • (vw) = efifi + en fa. 
Energy estimates on this equation then show that for all < T < T max , one has 

\ w \L°°([0,T/e],H<>) < l^C(T, \U\x%,)(\fl\L°°([0,T/e],H<>) + |/2|i«([0,T/e],JT*)) ■ 
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Now, one deduces from the explicit expression of fj (J = l,2)that \ fj\L°°([o,T/e],H<>) < 
C(\U\ x s+4.); with Theorem[TJ we deduce that \f 3 \l°°([o,t/£],h b ) — \Uq\ x s+4), 
and the result follows. □ 

Acknowledgments. The author is grateful to David Lannes for encouragement 
and many helpful discussions. 
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